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ABSTRACT 


Operational EW test and evaluation experiments require that the position of the 
aircraft and other moving objects on the range be known precisely as a function of time. 
Terminal Time-Space-Position Information (TSPI) systems involve the range platforms 
interacting at close distances and therefore require precise trajectory information over a 
restricted volume of space. Terminal TPSI systems are used for tactics evaluation and the 
evaluation of simulated weapons firings (e.g., captive-carry hardware-in-the-loop missile 
simulators). Distributed sensor TSPI systems consist of two or more measurement 
sensors located some distance from each other. Each sensor makes a measurement of 
target angle and range. Distributed sensor systems are more complex than single-point 
systems involving multiple hardware installations, complex mathematical computations 
to extract coordinate information, synchronization of multiple measurements and 
calibration of a number of different stations. 

This paper presents a novel distributed sensor TSPI architecture that provides 
precise positioning information of the target relative to a fixed inertial coordinate system. 
The architecture efficiently integrates the information from an inertial navigation system 
(INS), a global positioning system (GPS) and any number of distributed RF sensors 
which may be located onboard a captive-carry aircraft. The significance of this work is 
that by knowing the target's position in a fixed inertial frame of reference (derived from 
the integration process) an evaluation can be made as to the effectiveness of any 
electronic attack or off-board decoys that might have been launched during the field test 
scenario. The induced INS, GPS and sensor noise and the corresponding errors due to the 
integration process are evaluated numerically as a function of the weapon system being 
used. The accuracy in the targeting information is also quantified and compared with the 


true expected values. 
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I. INTRODUCTION 


A. CAPTIVE CARRY ASM EXPERIMENTS 


Currently, the Naval Research Laboratory, Washington, DC conducts captive-carry 
anti-ship missile (ASM) experiments onboard a modified P-3 aircraft. The purpose of the 
experiments is to evaluate the effectiveness of ship defensive systems against various 
types of anti-ship missiles. Currently, the evaluation process involves collecting 
information pertaining to the time and location of certain events which occur during a 
particular test scenario for post-test analysis. One of the goals of the post-test analysis is 
to accurately locate the P-3, the unit(s) involved in the exercise, and any offboard 
countermeasure devices (i.e., CHAFF, RUBBER DUCKS, etc.). By accurately locating 
the players in the scenario, the effectiveness of the countermeasures can be accurately 
determined. 

NRL's P-3 research aircraft is shown in Figure 1.1. It has the capability of carrying 
eight ASM Simulators beneath its main wings. Each ASM Simulator has an active seeker 
installed. Each ASM Simulator has the capability of providing the following information: 


e Target tracking accuracy 
e Seeker guidance type 

e Transmit frequency 

e Pulse width 

e Pulse repetition interval 
e Scan rate 


e Search sector 


e Beamwidth 
e Transmit Power 


e Antenna Gain 


e Azimuth, Elevation and Range to target 
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Figure 1.1: NRL's P-3 RESEARCH AIRCRAFT 


A typical scenario involves the P-3 aircraft equipped with the ASM Simulators and 
the appropriate data recording equipment flying towards a Naval Surface combatant with a 
self protection suite onboard (e.g., AN/SLQ-32(V)4). The P-3 makes several passes at the 
participating unit from several different aspect angles. Usually, the P-3 begins each run at 
a range of twelve nautical miles from the unit. At this range the seekers are turned on to 
begin the search for the target. As the aircraft approaches the ship, eventually, the seeker 


switches from search to track mode. In track mode the seeker locks on to the target in 


range and azimuth until the P-3 flies past the target or a countermeasure is employed. 


When a countermeasure is employed, the seeker may or may not break lock on the target. 


B. PRINCIPAL CONTRIBUTIONS 


The algorithms presented in this thesis provide a precise time-space-position 
information system (TSPI) based on the inputs from a number of distributed RF sensors, a 
GPS and an INS. The target's position in a fixed inertial frame of reference is derived 
from the integration of the INS/GPS/Radar systems that are located onboard NRL's P-3 
research aircraft. The significance of this thesis is that it provides the necessary tools to 
help NRL perform an evaluation of the shipboard self protection system effectiveness. By 
knowing the target's position as derived from the integration process an evaluation is also 
made as to the effectiveness of any launched decoys during the test scenario. These 
algorithms help the test engineers determine if the seeker broke lock on the target to track 
an offboard decoy (e.g. chaff). 

Also, the relevance of this thesis applies to other applications as well. By knowing 
a target's position relative to a fixed inertial reference frame, the targeting information can 
be sent to any platform and be understood immediately. This type of targeting information 
is ideal for remote firing of weapons that use GPS as an input. 

This thesis develops new algorithms to integrate the INS and GPS along with the 
distributed onboard seekers for determining the target position relative to a fixed inertial 
frame of reference. Mathematical models of the P-3, INS, GPS, and the ASM seekers are 
developed using SIMULINK. These models are used to provide the necessary inputs for 


the distributed sensor TSPI (DSTSPI) algorithms. Evaluation of the algorithms indicates a 


significant increase in the accuracy of the target position. The results make evident the 
importance of the INS/GPS integration process in knowing the ASM seeker's position in a 
fixed inertial frame of reference. The induced INS, GPS and RF sensor noise and the 
corresponding errors due to the integration process are evaluated numerically as a function 
of the weapon system being used. The accuracy in the targeting information is also 


quantified and compared with the true expected values. 
C. THESIS ORGANIZATION 


Parts of this thesis assume the reader has a background in the areas of radar and 
aeronautics. Thorough discussions are provided in the areas specifically relevant to the 
algorithms presented in this thesis. 

In Chapter II, the background information needed to gain a full appreciation of this 
thesis is discussed. First, the errors associated with the Global Positioning System (GPS) 
and Inertial Navigation System (INS) are discussed along with basic system descriptions. 
Because both GPS and INS reside within fixed inertial frames of reference, a thorough 
discussion of the fixed inertial coordinate systems that are used in this thesis are provided. 
Also, because of the dynamics of the P-3 during the test scenario, a general methodology 
on how to use lever arm corrections is discussed. Finally, this chapter discusses the 
equations used in the modeling of the monopulse and conical scan seekers that are used in 
the simulation. 

In Chapter III, the DSTSPI algorithms are thoroughly developed. Here, a general 
description of the process is given along with an explanation of the SIMULINK block 


diagrams that are used. 


In Chapter IV, a discussion is provided on the results of the DSTSPI algorithms. 
Several plots are used to aid in the validation and analysis of the results. In addition, a 
further examination 1s made as to the utility of the algorithm and whether or not it can live 
up to the user’s expectations based on the relevance to the ASM/EW test scenario. 

Chapter V discusses some important considerations involving the hardware and 
software integration process onboard the P-3. Finally, Chapter VI discusses conclusions 


that are made in the process of this analysis. 
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Il. BACKROUND INFORMATION 


A. GLOBAL POSITIONING SYSTEM 


The Global Positioning System (GPS) consists of 24 satellites that orbit the earth. 
Of the 24 satellites three are ready service spares. GPS provides locating information to 
an unlimited number of users. Today, GPS is widely used by both the military and the 
civilian sector. More and more civilians are purchasing inexpensive GPS receivers for 
recreational use. The applications of GPS seem limitless. [Ref. 1, p. 121] 

GPS is no doubt the most powerful navigational aid available today. There are, 
however, a few drawbacks to the system. First, the positioning information available from 
GPS is updated at a 10 Hz rate for high-end systems and at a 1 Hz rate for most other 
systems. The update rate becomes important when the application involves real time 
processing. In the case where an update rate greater than 10 Hz is required, the 
positioning information from the GPS receiver must be smoothed by employing a Kalman 
filter or some other predictive algorithm. Some higher-end GPS receivers do provide a 
smoothed output that the user is able to sample at higher frequencies. 

There are several factors that determine the GPS accuracy of each position fix. 
The accuracy can vary based on the time of day, the location of the GPS antenna relative 


to the GPS satellites and the type of GPS receiver used. 


The major sources of error for GPS are summarized as follows: [Ref. 1,p 125] 


e Selective availability 
e Clock differences 

e Receiver noise 

e Multipath 

e Atmospheric delays 


e Ephemeris error 

The sources of error are detailed below. 

1. Atmospheric Delay 

Atmospheric delays occur mainly in the ionosphere and the troposphere. In the 
ionosphere, time delays occur because the charged particles interact with the transmitted 
GPS signals. As a result of this interaction, the GPS signal arrives at the receiver with 
some inherent time delay. Since the rate of propagation is nearly the speed of light, any 
delay in the receipt of the signal translates to an increase in the calculated range. 

The equation for the ionospheric delay is given by 


ie 40.3 
Siar 





TEC (2-1) 


where At is the delay in seconds, f is the frequency in Hz, ¢ 1s the speed of light in m/s, 
TEC is the total electron content (electrons/m’) along the signal path and 40.3 is an 
empirically derived constant with units (m’/s’/electrons). Since, in the ionosphere, the 
charged particles are typically between 100 to 1000 kilometers above the earth, the path 
length for interaction can be very large. The pseudorange error can be as great as 40 
meters depending upon the time of day, time of year, the solar cycle and the geographic 


location. [Ref. 1,p 126] 


The lower part of the atmosphere can also cause significant delays in the GPS 
signal. This portion of the atmosphere is called the troposphere. These delays are caused 
by the interaction of the GPS signal with water molecules entrained in the atmosphere. 
Here, a 2 meter error can occur in the pseudorange calculation when a satellite is directly 
overhead and a 28 meter error when a satellite is just above the horizon. There are models 
available that can predict the errors from both the ionospheric and tropospheric delays. 
These models improve the performance of the GPS receiver. 

2. Selective Availability 

Selective availability is the status of GPS when the Department of Defense 
degrades the position accuracy to 100 meters. This is accomplished by dithering the 
individual satellite’s clock signal. Dithering is a process by which the clock signals of the 
satellites are encoded such that only military users are able to take advantage of the 
highest accuracy. 

3. Clock Errors 

No matter how precise and accurate a clock is said to be, all clocks are plagued by 
bias and drift errors. Bias is a fixed error associated with the offset of the clock from the 
correct time. Drift is defined as the rate of change of the accuracy of a clock with time. 
Typical clock errors are on the order of a few nanoseconds. A one nanosecond error 


equates to roughly three tenths of a meter error in the pseudorange calculation. 


4. Ephemeris Error 

Ephemeris errors occur when a satellite’s position is different from the position 
contained within the navigation message that is sent by each satellite. The Ephemeris 
errors are typically on the order of a half meter in the radial direction and three and a half 
meters in the planar direction. 

5. Multipath 

Multipath errors result when the signal does not take a direct path from the satellite 
to the receiver. A signal may bounce off a particular object before reaching the GPS 
antenna. This causes delays in the reception of the signal which can directly translate into 
errors in the pseudorange calculation. Because of the design of current GPS receivers, this 
effect is usually insignificant and is ignored. 

6. Receiver Noise 

All receivers inherently add noise to the signal they receive. A typical receiver has 
a noise value which contributes around 7.5 meters to the error in the pseudorange 
calculation. The overall error associated with receiver noise is wholly dependent upon the 
quality of the receiver being used. 

One method in eliminating most of the errors inherent with GPS is to use 
Differential GPS. Differential GPS (DGPS) involves having another GPS receiver fixed 
to a known location within close proximity of the user’s GPS receiver. This is called a 
differential GPS station. This differential station becomes susceptible to the same errors 
as the user’s GPS receiver. Because the exact location of the differential station is known, 


the differential station can transmit its location along with the corrections for the errors 
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discussed above. This process greatly improves the accuracy of the positioning system. 
An advertised value for position accuracy for a DGPS receiver is around one meter. 
[Ref. 3] 

The DSTSPI algorithms use a DGPS model for the input of the GPS positioning 
information. This model assumes a position accuracy of one meter with an update rate of 
10 Hz. The development of this model is not discussed in this work. Details are found in 


Reference 1, pages 141-152. 
B. INERTIAL NAVIGATION SYSTEM 


The inertial navigation system (INS) detects changes in the total acceleration of a 
body relative to an inertial frame of reference. The primary component of the INS is the 
inertial measuring unit (IMU). The IMU consists of three accelerometers, three rate gyros, 
and two inclinometers. The accelerometers measure the linear, centripetal and 
gravitational acceleration effects on the body. The INS removes the gravitation 
component of the acceleration from the total acceleration. This resultant acceleration is 
integrated once over time to find the velocity of the body and twice to find the position of 
the body. The rate gyros of the IMU measure the angular velocities of the body. The 
inclinometers measure the orientation of the body relative to the inertial frame of 
reference. This orientation is typically denoted as the pitch, roll and yaw. [Ref. 1, p. 136- 
Pe) 

There are two types of inertial navigation systems. The first type is the gimbaled 
system. Here, the gimbaled IMU is allowed to rotate about four gimbals while a controller 


maintains the IMU platform in a local level orientation towards true North. The relative 


differences between the gimbals and the platform provide the inertial quantities. The 
problem with the gimbaled systems is that they are usually very heavy and bulky. 
Because of this, the gimbaled system is not very well suited for an aircraft 
implementation. 

The second type of inertial navigation is the strapdown system. Here, the IMU is 
literally strapped down to the aircraft. In this manner the IMU maintains a constant 
orientation to the local coordinate system of the aircraft. Because of this configuration, the 
output of the IMU is given relative to the aircraft’s coordinates. The INS transforms the 
output of the strapdown IMU to accelerations and angular velocities relative to the inertial 
frame. These transformations require substantially more computing power than what is 
needed for the gimbaled system. Because computers are relatively fast and inexpensive, 
the strapdown system is the one generally found on most aircraft. 

The error sources for the INS are mainly due to errors pertaining to the IMU. 
These errors are categorized as follows: bias, cross-axis sensitivity and noise floor. The 
bias error results when the output has a constant offset from the correct value. If there are 
no accelerations on the body, the total acceleration should be zero. Any other value would 
indicate a bias error. A cross-axis sensitivity error occurs when the IMU 1s not properly 
aligned with the aircraft’s coordinate system. Because it is virtually impossible to achieve 
a perfect alignment of two mechanical systems, a cross-sensitivity error is always present. 
The noise floor for INS prevents measurements below a specified value. It acts as a 
threshold for the output of the INS. The noise floor is caused by the mechanical and 


electrical components of the INS. 


The DSTSPI algorithms use a strapdown INS model for the input of the pitch, roll 
and yaw of the P-3. This model assumes an accuracy of 0.5 degrees in the determination 
of the pitch, roll and yaw. The development of the INS model is not discussed in this 


work. Details are found in Reference 1, pages 152-157. 


C. COORDINATE SYSTEMS 


Several coordinate systems are used in the development of the DSTSPI algorithms. 
For this reason, a brief discussion is provided on the basic characteristics of these 


coordinate systems. 
1. Earth-Centered-Earth-Fixed 


The Earth-Centered-Earth-Fixed (ECEF) coordinate system is a mght-hand 
orthogonal coordinate system that is fixed at the Earth’s origin [0, 0, 0]. Figure 2.1 shows 


the ECEF coordinate system in relation to the Earth. 


Z 


PRIME MERIDIAN 





Figure 2.1: Earth-Centered-Earth Fixed Coordinates 
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The x-axis is defined as the direction from the origin through the intersection of 0° 
longitude and the equator. The y-axis is defined as the direction from the origin through 
the intersection of 90° east longitude and the equator. The z-axis is defined as the 
direction from the origin through the point defined as true north, or 90° latitude. The 
ECEF coordinate system rotates with the earth. One important advantage of the ECEF 
coordinate system is that all positions that are referenced to ECEF are independent of the 


model used for the irregularities of the Earth’s surface. [Ref. 2, p. 4] 
2. Geodetic Coordinate System 


The geodetic coordinate system is the coordinate system that is most commonly 
used for navigation. Generally, objects which are fixed to the Earth’s surface are located 
in the geodetic coordinate system by designating their latitude, @, and longitude, A. The 
latitude and longitude are defined by their relation to the prime meridian and the equator. 
The geodetic latitude, 0, is defined as the angle from the equatorial plane to the vertical 
direction of a line normal to the reference ellipsoid passing through a location on the 
reference ellipsoid. The geodetic longitude, A, is defined as the angle between the 
reference plane and a plane passing though a location on the reference ellipsoid defined by 
N, where both planes are perpendicular to equatorial plane. The geodetic height, h, is 
defined as the distance from the reference ellipsoid to the location in a direction normal to 
the reference ellipsoid. The reference ellipsoid is a model that defines the shape of the 
Earth. Most modern navigation systems use the World Geodetic System -- 84 model 
(WGS-84). The WGS-84 model gives the model parameters for the shape of the Earth 


along with the local height correction for a particular latitude and longitude. [Ref. 2, p. 5] 
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Figure 2.2 shows a general depiction of the geodetic coordinate system. Here, N 
is the distance from the polar axis to the surface of the reference ellipsoid, h is the height 


correction for the terrain on the Earth’s surface and P defines a position in geodetic 


coordinates. 
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Figure 2.2: Geodetic Coordinate System 


3. Tangent Plane Coordinate System 


The tangent plane coordinate system is a nght-hand orthogonal system. The origin 
of the tangent plane coordinate system is defined by passing a plane through a point on the 
earth’s surface which is defined by a particular latitude and longitude. The axes, (X,Y,Z), 
of the tangent plane coordinate system are defined as: [Ref. 2, p. 5] 


e X-axis points toward true North 
e Y-axis points toward true East 
e Z-axis points in the downward direction from the point where the tangent plane 


intersects the reference ellipsoid to a point on the polar axis on the earth 


Figure 2.3 shows a general depiction of the tangent plane coordinate system. 





Figure 2.3: Tangent Plane Coordinate System 


4. Navigation Coordinate System 


The navigation coordinate system is a right-hand orthogonal coordinate system 
whose origin is collocated with the aircraft’s center of gravity. It maintains a local level 
orientation to the reference ellipsoid in the same manner as the tangent plane. If the origin 
of the aircraft is collocated at the ongin of the tangent plane, then the navigation 
coordinate system and tangent plane coordinate system are the same (assuming the 
orientations of the axes are the same.) The origin of the navigation coordinate system 1s 
defined to be the center of the aircraft’s inertial navigation system. The orientation can be 
chosen to be East North Up (ENU), NED, etc.. [Ref. 2, p. 6] 

Figure 2.4 shows the navigation coordinate system with a North East Down (NED) 
orientation. This is the orientation used in the DSTSPI algorithms developed in this 


thesis. It 1s convenient to choose NED for the onentation of the axes because the 


coordinate system naturally aligns itself with true North. The difficulty in dealing with 


NED orientation 1s that altitude is a negative quantity. 





Figure 2.4: Navigation Coordinate System 


5. Body Fixed Coordinates 


The body fixed coordinate system is the nght-hand orthogonal coordinate system 
that is attached to the aircraft and the axes of which are defined by the input axes of the 
inertial navigation system. The body fixed coordinate system rotates with the aircraft 
which implies that the body accelerations of the aircraft are the same as the body fixed 
accelerations if the inertial navigation system is a strapdown system. Figure 2.5 shows a 
depiction of the body fixed coordinate system on the P-3. The orientation of the axes are 
forward, right wing and down. The orientation is defined as a matter of convention. 
Notice that this orientation is the same orientation used for the navigation and tangent 
plane coordinate systems (e.g., forward in the body fixed coordinates is North in the 
tangent plane coordinates.) Keeping the axes orientations the same for these coordinate 


systems results in simpler coordinate transformations. [Ref. 2, p. 8] 
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Figure 2.5: Body Fixed Coordinate System 


D. LEVER ARM CORRECTION 


To accurately locate the seekers for each ASM Simulator relative to a fixed 
coordinate system such as ECEF, the first problem that must be solved is the fact that each 
seeker is not collocated with the GPS antenna. If they were, the seeker’s location would 
simply be the output of the GPS receiver. Since positioning information from the GPS 
receiver 1s given relative to the GPS antenna in ECEF or geodetic coordinates, natural 
position errors are inherent based upon the displacement between seeker and GPS antenna. 
In addition, because the P-3 is not, in general, locally level with the earth’s surface (the 
reference ellipsoid), additional errors are inherent because of the pitch, roll and yaw of the 
aircraft. To correct for the position errors that result from GPS antenna displacement from 
each seeker and from the pitch, roll and yaw of the aircraft, a lever arm correction m 1s 


used. [Ref. 3, discussion] 


The general algorithm for performing a lever arm correction is shown in Figure 
2.6. Three things must be known before a lever arm correction is made. The first is some 
known fixed inertial position that is located somewhere on the P-3. This is usually the 
position of the GPS antenna given in ECEF coordinates, or the position of the INS given 
in navigation coordinates. The fixed inertial coordinates that are chosen are those in 
which the intended local body position is to be referenced (e.g., if the user wants to 
reference all positions to the tangent plane coordinate system, then this is the coordinate 
system that should be chosen). This position is denoted ™P,,, where the subscript, PT, is 
defined as the local position (local positions refer to those positions located in the body 
fixed coordinate system which are to lever arm corrected to a fixed inertial coordinate 
system.) An example of a local position is the seeker location. The superscript, FI, is 


defined as the fixed inertial coordinate system in which the local position 1s to be resolved. 
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Figure 2.6: Lever Arm Correction Algorithm 
To accurately perform a lever arm correction, the positions of the GPS antenna and 


the individual seekers for each ASM Simulator must be known relative to the center of the 


aircraft in the body fixed coordinate system. These measurements are extracted from the 
existing drawings of the P-3 or taken manually. The overall accuracy of the lever arm 
correction depends upon the accuracy and precision of the measurements of the relative 
positions of the distributed sensors. Each simulator has two position vectors that must be 
determined in body fixed coordinates. One position vector describes the location of the 
associated distributed seeker and the other describes the after position of the ASM 
simulator. These two position vectors are used later in the analysis to determine the 
attitude of the Simulator. Notice that the forward location of the ASM simulator is 
essentially the location of the seeker. These position vectors are notionally depicted in 


Figure 2.7. 
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Figure 2.7: General layout of systems requiring lever arm corrections 


Once the position vectors for the GPS antenna and the ASM Simulators are known 
in the body fixed coordinate system relative to the origin of the P-3 (which is the location 
of the INS), the next step is to compute the rotation matrix between the body fixed 
coordinate system and the tangent plane coordinate system’. This is done by first knowing 


the pitch, roll and yaw of the aircraft. The pitch, roll and yaw of the aircraft are taken 
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from the inertial navigation system (INS). The pitch, roll and yaw are defined as the Euler 
angles between the body fixed coordinate system and the tangent plane coordinate system 
when the origins of the two coordinate systems are collocated. We define pitch as ©, yaw 
as ‘Y, and roll as ®. [Ref. 1, p. 17] 

In order to derive the rotation matrix between the body fixed coordinate system 
and the tangent plane coordinate system the following notation and assumptions are used: 


e {T} represents the inertial tangent plane coordinate system. 

e {B} represents the body fixed coordinate system. 

e The INS is located at the center of gravity of the P-3 aircraft or the origin of 
{B}. Any displacement errors of INS from the true center gravity of the P-3 are 
assumed to be negligible. 

e The {B} coordinate system is defined with X, as the direction forward of 
motion with Y, pointing in the starboard direction (right wing) and Z, pointing 
in the downward direction. 

e The {T} coordinate system is defined with X, in the direction of true north, Y, 
in the direction of east and Z, pointing down towards the center of the earth 

e Roll, ®, is positive when the rotation about X, 1s clockwise when looking in 
the positive X, direction. 

e Pitch, ©, is positive when the rotation about Y, 1s clockwise when looking in 
the positive Y, direction. 

e Yaw, ‘Y, is positive when the rotation about Z, is clockwise when looking in 
the positive Z, direction. 


The Euler angles are depicted in Figure 2.8. Using the nght-hand rule defines the 


direction for rotation for roll, pitch and yaw. 


"Note that the tangent plane and navigation coordinates systems differ only by the altitude of the aircraft. 
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Figure 2.8: Directions Euler angles relative to body fixed coordinates 


When finding the rotation matrix, the order in which the rotation is taken is 
important because Euler angles are not vectors. The sum of ®+@ is not generally equal to 
©+@. For purposes of this derivation the order is taken to be Z-Y-X, or rotate about Z 
first, then Y and X last. [Ref. 1, p. 18] [Ref. 2, p.13] 

It is convenient to collocate the origins of the body fixed and tangent plane 
coordinate systems as seen in Figure 2.9. By collocating these coordinate systems, the 


trigonometric relationships between them are easily found. 
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Figure 2.9: Collocated body fixed and tangent plane coordinate systems 


Given a position vector, "P, which is resolved in {B}, the goal is to find the components of 
5P that are resolved in {T}. To do this, first rotate the {B} coordinates about Z,. The 
angle by which {B} must be rotated is defined by ‘¥. Therefore, the first rotation matrix is 


a function of ‘¥. Figure 2.10 shows the rotation about the Z, axis. 
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Figure 2.10: {B} rotated about Z, 
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After the rotation through the yaw angle the new axes are defined as ie Ves Z',| 


The tngonometric relationships between the unprimed and primed frames are used to 


derive the transformation matrix between the two frames. These trigonometric 


relationships are depicted in Figure 2.11. 
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Figure 2.11: Trigonometric relationships when rotated about Z, 


From Figure 2.11, we see that the components of the new primed system are a 
function of the unpnmed system and ‘VY. X"', is simply the result of the projection of X, 
onto X', and the Y,sin'¥ term. This results in Equation (2-2). Equation (2-3) is derived in 
a similar fashion. Also, since the rotation occurred about the Z, axis the z components do 
not change during the rotation. This is summarized by the following relationships: 

X, = X, cos Y + Y, sin (2-2) 
Y, =—X, sin VY + Y, cos V (2-3) 
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Ae (2-4) 


In matrix form equations (2-2), (2-3) and (2-4) simplify to 


Xz} | cos¥ sin’ 0] x, 
Y, |=|—sin¥ cos¥ 0] Y, (2-5) 
7 0 0 1 | Z, 
Now define the rotation matrix for ‘¥ as 
cos ¥ sin? 0 
T(Y)=|-sin¥ cos¥ 0]. (2-6) 
0 0 1 


Next, rotate the primed coordinates about Y’, through the angle ©, which is the 


pitch angle of the aircraft. This rotation is seen in Figure 2.12. 
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Figure 2.12: {B} rotated about Y’, 


The trigonometric relationships for this rotation are defined 1n Figure 2.13. 
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Figure 2.13: Trigonometric relationships when rotated about Y’, 


From Figure 2.13, the double primed coordinates is defined by the following 


relationships: 
X, = X, cosO —- Z, sinO (2-7) 
eae (2-8) 
Z,, = X,sin©+Z, cos© (2-9) 


In matrix form these equations become 


X,| [cos®© 0 -sinO]| x, 
Y, |=| Ome" 1 0 ae (2-10) 
Z., sin©® 0 cos© | Z, 


with the rotation matrix for © as 
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cos® 0 —-sinO 
T(O)=| 0 1 0 |. (Cain); 
sin® 0 cos® 


Next, rotate the double primed coordinates through the angle ®, which is the roll 


of the aircraft. The result is depicted in Figure 2.14. 
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Figure 2.14: {B} rotated about X”’, 


The trigonometric relationships between the double primed and triple primed coordinates 


are depicted in Figure 2.15. 


20 


y"', cos 


: +X Pid * 
YY", v OUT OF PAGE 





* 


- ¥", sin ® 


Figure 2.15: Trigonometric relationships when rotated about X”, 


From Figure 2.15, the triple primed coordinates are defined in terms of the double primed 


coordinate and the angle ®. The result is the following equations: 


ee (2-12) 
Y, = Y, cos + Z, sin® (2-13) 
Z,, =~Y, sin® + Z, cos® (2-14) 


In matrix form, these equations become 


X;| [1 0 0 |x, 
Y, |=|9 cos®  sin® jj Y, |, (2-15) 
Z. 0 ~—sin® cos® || Z, 


with the corresponding rotation matrix for ® defined as 
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1 0 0 
T(®)=|0 cos® sin®]. (2-16) 
0 -sin® sin® 
Combining Equations (2-5), (2-10) and (2-15) and adhering to a [Z-Y-X] rotation 


sequence results in the following equation: 


xX 1 0 0 cos® Q —sin®| cos¥ sin? O}] X 
Y!| =/0 cos® = sin® 0 1 0 —sin¥ cos¥ OY; (2-17) 
Z Pi 0 -—sin® cos®M]|sin©® 0 cos® 0 0 1} Z 


T) 

For a lever arm correction, the lever arm coordinates are known in the body fixed 
coordinates. Equation (2-17) applies to knowing a position vector in tangent plane 
coordinates and then finding the position 1n the body fixed coordinates based on the Euler 
angles. To find the lever arm position vector resolved in the tangent plane coordinates, 


find the inverse of Equation (2-17). The matrix form that results is conveniently written as 


xX 1 0 0 cos® 0 —sin®} cos? sin O]}}]) X 
Y| =INV4/0 cos® sind} 0 1 #40 |—sin® cos¥ Ol} Y| . (2-18) 
Z 0 -sin® cosD]sn© Q = cos® 1 0 1|}| Z 
{T) {B) 
Now, define the rotation matrix, ,R as 
1 0 0 cos©® Q -—-sin®| cos¥ sin? 0 
;R=INVs|0 cos® sin®| 0 1 0 —sin¥ cos¥ O|7. (2-19) 
0 -sin® cosM]|sin® QO = cos® 1 0 1 


which 1s the rotation from the body fixed coordinates to the tangent plane coordinates. 


This 1s more compactly written as 


"R = INV{T(®)T(O)T(Y)} . (2-20) 
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Now, by taking a known position that is resolved in tangent plane coordinates, say 
the GPS antenna position, and adding this position to a position vector that is rotated from 
the body coordinates to the tangent plane coordinate, the result is a position in the body 
fixed coordinate system that 1s resolved in a fixed inertial frame of reference. To find a 
position resolved in tangent plane coordinates given a position resolved in body fixed 
coordinates with the origins collocated the following equation is applied 

P= hoe (2-21) 
where the subscript 1 means a point somewhere in the body fixed coordinate system. 
Now, by adding this position to the GPS position in the tangent plane coordinate system 
(fixed in the body frame), results in 

*P = 'PopstpR*P,. (2-22) 
which is the addition of two position vectors resolved in the same coordinate system. 


Equation (2-22) is a simple recipe for performing a lever arm correction. 
E. SEEKER MODELS 


This thesis uses both a conical scan on receive only seeker (COSRO) model and a 
monopulse seeker model to help in the validation process of the DSTSPI algorithms. 
These models are used because they reflect the most common types of seekers found in a 
typical missile inventory throughout the world. The parameters chosen for these models 
are based on typical performance parameters associated with each seeker type. These 


parameters are summarized in Table 2.1: 
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Estimated Parameter 
Peak Power, P; (kilowatts) 
First side lobe Ratio (dB) 
Frequency, f (GHz) 
Antenna Gain, G (dB) 


1 
2 


Noise Bandwidth, B (MHz) 
Noise Figure, F (dB) 


0 


Range resolution (meters) 





a 
= 
Half Power Beam Width, 0,,, (degrees) = 


Number of Pulses Integrated, n 


Table 2.1: Estimated seeker parameters 


1. Single Pulse Signal to Noise Ratio 


The first step in modeling a seeker 1s to estimate the signal to noise ratio (S/N) as a 
function of the range between the seeker and the target. First, an assumption must be 
made about the radar cross section of the target, o. In general, the radar cross section of a 
typical naval surface combatant is of the order of 3000 m’. [Ref. 5, p. 129] 

The peak power radiated by any radar transmitter is defined as P, in watts. The 
power density at a given range, R, from an omnidirectional antenna is given by [Ref. 7] 


and [Ref. 4, p. 10] as 


Pe“ 
4nR’ 





(2223) 


where « is the attenuation coefficient given in nepers/km. The attenuation coefficient is 


estimated to be 0.055 dB/km at both 9 and 17 GHz. [Ref. 4, p. 278] 


3] 


The power density from a directional antenna with an antenna gain of G, is given by 


Ge — 


4nR’ a) 


The power that is intercepted by the target and re-radiated back to the radar antenna is 


given by 
P,G,ce ™* 
——— 2-25 
4nR’ ( ) 
where o 1s the radar cross section of the target in square meters. [Ref. 4, p. 11] 
The power density received from the re-radiated power is given by 
P.G.oe”™ 
tt — (2-26) 
(4nR 
since the attenuation and spreading are over twice the distance. 
The signal collected by the antenna is given by 
P.G.A,oe* 
=e = (2-27) 
(4x) R 
where A, is the effective area of the receiving antenna which is equivalent to 
Gar 
Ae (2-28) 
4n 


Now, by assuming the transmitted gain of the antenna, G,, 1s equivalent to the received 


gain of the antenna, G,, the received signal power 1s approximated by the following: 


2 ~2aR 4 2 
P,G,oe "A 


(ane (watts). (2-29) 


Now, define the noise power in the receiver as 


Sy 


NOISE = kTBF , (watts), (2-30) 
where k is Boltzman's constant, T 1s the temperature 1n degrees Kelvin, B is the receiver 
noise bandwidth and F is the receiver noise figure. The temperature is assumed to be 300 
degrees Kelvin. Therefore, kT is 4.14 x 10” watts/Hz at 300 K. The single pulse signal 
to noise ratio is defined as the signal power collected by the antenna divided by the noise 
in the receiver or 


P,G,oe "2? 


= any a eB) 


(22200) 


Given the single pulse signal to noise ratio as a function of range, the precision in 


the angular measurement for a COSRO and monopulse seeker can be evaluated. 

2. COSRO Seeker 

For a COSRO type seeker, the error caused by thermal noise over the n pulses 
integrated by the servo loop in each coordinate (azimuth and elevation) is [Ref. 4, p. 383] 


Osan L, 


= 280+ 232 
ie k,(S/N)_n ne 


where k, is the conical-scan error slope, (S/N),, 1s the single pulse signal to noise ratio, n 


is the number of pulses being integrated and L, is the crossover loss due to beam squint. 
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Figure 2.16: One-way error slope [Ref. 4, p.389] 


The conical scan error slope, k,, and the crossover loss, L,, are shown in Figure 
2.16 as a function of the normalized squint (offset) angle, 0,/03,,. It 1s assumed that the 


beam has a Gaussian shape which results in an optimum error slope of 1.0 at a normalized 


offset angle of 0.5. 


With &, / L, =1, the error due to thermal noise becomes 


Os58 
0, = oS. 2-33 
: J(S/N)_n es 


The scintillation error associated with a conical scan seeker 1s given by 


) 
g.= 022s{ Sx " (2-34) 
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where t, is the envelope correlation time, f, is the nutation frequency, and B, is the servo 
bandwidth. Assuming a target with an envelope correlation time t..= 0.06 s, a scan 
frequency of f, = 60 Hz and a servo bandwidth 8, = 2 Hz, for an optimized squint angle 
0, = 0.50,,, results in a scintillation error of 
o, = 0.01250,,,. (2-35) 

The total error in the precision of the angular measurement for a conical scan 
seeker 1s found by adding Equations (2-33) and (2-35). The error in the angular 
measurement, 60, or 60 is assumed to be Gaussian distributed and is given by 

1( uve) 

_ 


oe 
2n(o? +05) , 


(2-36) 
where o, 1s a function of the single pulse signal to noise ratio. These equations and the 
parameter estimations found in Table 2.1 are used to develop a model for a COSRO seeker 
type. 


3. Monopulse Seeker 


The next step is to find the precision of the angular measurement for a monopulse 
seeker. The precision of the angular measurement in a thermal noise environment is 


0 3dB 


Somp > Kk J2(S/N)n (2-37) 


Here, k,, 1s the difference slope and n is the number of pulses being integrated. The value 


for the monopulse error slope is determined from Figure 2.17. It is assumed that the shape 
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of the pulse is Gaussian. From Figure 2.17 with a first sidelobe of 24 dB, the value of k,, 


is approximately 1.9. 
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Figure 2.17: Difference slope versus sidelobe level [Ref. 4,p. 403] 


The error in the angular measurement, 50, or 6d is assumed to be Gaussian 


distributed and is given by 


(2-38) 
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where ©, is a function of the single pulse signal to noise ratio. These equations and the 
parameter estimations found in Table 2.1 are used to develop a model for a monopulse 


seeker type. 


oy 





Hil. ALGORITHM DEVELOPMENT 


A. GENERAL DESCIPTION OF ALGORITHM 


This section provides a general description of the DSTSPI algorithms. One of the 
goals of the DSTSPI algorithms is to minimize the errors in the derived target position due 
to the displacement of the GPS antenna from the positions of the distributed ASM seekers 
and the pitch, roll and yaw of the P-3. As mentioned in Chapter I, this is accomplished by 
using lever arm corrections. Another goal of the DSTSPI algorithms is to provide the 
necessary coordinate transformations and corrections to the seeker-to-target information 
resulting in a derived target position in a fixed inertial coordinate system. 

Figure 3.1 depicts the overall functionality for the algorithm. The integration 
algorithms perform the necessary lever arm corrections and compute the pitch and yaw of 
each Simulator. The output from the integration algorithms in conjunction with the target 
data received from the ASM seekers 1s used to derive the target position in a fixed inertial 


frame of reference. 
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Figure 3.1: General DSTSPI functionality 


B. EXPLANATION OF SIMULINK BLOCK DIAGRAMS 


In the development of the DSTSPI algorithm, SIMULINK® is used as the primary 
tool for the analysis. SIMULINK is a non-linear simulation package developed by The 
MathWorks Inc. for use with MATLAB®. One distinct advantage in using SIMULINK is 
that the results can be viewed in real time or sent to a file for post simulation analysis. 
Another distinct advantage of SIMULINK is that, like functional block diagrams, it 1s easy 
to see the logical progression of the algorithm. For this reason several SIMULINK block 


diagrams are included in this section. 
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The underlying assumptions for the algorithm are as follows: 


e GPS antenna is located at the [-S5m, 0m, 0m] position in the body fixed 
coordinates and its position is denoted by Paxr. 

e INS ts located at the center of gravity of the P-3 which is the origin of the body 
fixed coordinate system and its position is denoted by Pen. 

e Each Simulator is assumed to be 4 meters in length, which is strictly arbitrary. 

e Simulator #1 is located between the body fixed coordinates [2m, 5m, 0m] and 
[-2m, 5m, 0m]. The aft position of the Simulator #1 is denoted by P,,. The 
forward position of the Simulator or the seeker is denoted by P,,. 

e Simulator #2 is located between the body fixed coordinates [2m, -5m, 0m] and 
[-2m, -5m, 0m]. The aft position of the Simulator #2 is denoted by P,,. The 
forward position of the Simulator or the seeker is denoted by P,,. 

e The target position is denoted by Pyap. 

e The derived target position is denoted by P,,,,. 


e DGPS is used with an update frequency of 10 Hz. 


Figure 3.2 shows the relationships of the ASM Simulators and the GPS antenna in the 


body fixed coordinate system. Here, the values for the positions are given in meters. 
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Figure 3.2: Sensors in the body fixed coordinates (positions given in meters) 
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The notation used in the development of the DSTSPI algorithms is unique and 


needs to be well defined. The following is a list of notational definitions: 


e A superscript in front of the position signifies the coordinate system the 
position is resolved in. For example, *P,,7 is the GPS antenna resolved in the 
body fixed coordinates and 'P,,,7 is the GPS antenna position resolved in the 
tangent plane coordinate system. 


e The rotation matrix is denoted as ¥R where X is the coordinate system from 


which the rotation is occurring and Y is the coordinate system to which the 
rotation is to occur. 

e A primed variable indicates a value with no associated noise added to it. These 
are the actual or true values associated with the flight of the P-3 (e.g., position, 
Euler angles, etc.). 

e An unpmmed variable indicates a value with noise added to it. These are the 


derived values (e.g. GPS, INS, etc.). 
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Figure 3.3: DISTRIBUTED SENSOR TSPI Block Diagram 
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Figure 3.3 is the top-level block diagram for the DSTSPI algorithms. Here, the 


“Pp-3 FLIGHT PROFILER?” provides the GPS position at the aircraft’s center of gravity 
("Porn ), the INS pitch, roll and yaw(©,®,'¥ ), the true position of the aircraft (* Puen ) 
and the true pitch roll and yaw of the aircraft (© ,®,'Y ). Because the GPS position is 
the location of the GPS antenna in tangent plane coordinates, the block “CONDITION 
GPS” converts the center of gravity position ("P...) to a position located at the GPS 
antenna ('P,.,). The “TARGET MOTION?” block generates a target in tangent plane 


coordinates ("P_,.) for use in the algorithm and the validation process. The 
“SIMULATOR #1” block is identical to the “SIMULATOR #2” block except for the 
specific values used in the lever arm corrections and the seeker models employed. For this 
reason, the discussion of the DSTSPI algorithms only includes Simulator #1. The 
components of the “SIMULATOR #1” block are discussed in more detail later in this 


chapter. 
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Figure 3.4: P-3 FLIGHT PROFILER Block Diagram 


The “P-3 FLIGHT PROFILER” is depicted in Figure 3.4. The development of this 
part of the model was performed by LT Chuck Gill, et al, as part of a class project in the 
Aeronautical Engineering Department at the Naval Postgraduate School. A thorough 
discussion of this part of the algorithm is found in Reference 6. Although this block is an 
integral part of the analysis, its only function is to provide the outputs required for this 


algorithm which are assumed to be available on the P-3 test flight vehicle. 
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Figure 3.5: P-3 FLIGHT PROFILER CONTROLLER Block Diagram 


The controller used with the P-3 flight profiler required a few modifications to 
make it fully compatible with the DSTSPI algorithm. First, the controller, which is shown 
in Figure 3.5, 1s modified to allow the aircraft to steer courses in both the starboard and 
port directions. Prior to this modification, the aircraft would only tum in the port direction 
regardless of the ordered course. This correction is made by writing two MATLAB 
programs which corrects the inconsistencies in the flight controller. These programs are 
listed in Appendix A as “HEAD.M” and “OVER2PI.M.” Another problem corrected is 
the instabilities in the controller for steer commands greater than, say, 30 degrees. This 


correction is made by limiting the difference between the ordered steer course and the 
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actual steer course. This gave the controller time to adjust its output to match to required 


response. 
SF 
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Figure 3.6: TARGET MOTION Block Diagram 


The “TARGET MOTION?” block diagram is shown in Figure 3.6. The target 
motion is simulated by assuming an initial displacement from the location of the origin of 
the model’s tangent plane coordinate system (which can be any fixed inertial coordinate 
system). The initial displacement of the target is [5, 0, 0] in nautical miles, with the x 
axis pointing toward true North, the y axis pointing towards true East and the z axis 
pointing downward. The velocity of the target is fixed at 22 knots in a westerly direction 
which translates to -37.131 feet/second in the y direction. By simply multiplying this 
fixed rate by the elapsed time a new target position is calculated for each time step. The 
true position of the target ("P,,,) is then stored in the file “CST.MAT” for post 
simulation analysis. The block is changed to simulate different types of target motion. 


This is done to obtain the results in Chapter IV. 
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Figure 3.7: CONDITION GPS Block Diagram 


The “CONDITION GPS” block diagram shown in Figure 3.7 performs a lever 
arm correction to translate the GPS position at the center of gravity of the aircraft (* Popy ) 
to the GPS position at the GPS antenna location (*P,,,). As mentioned earlier, the GPS 


antenna is assumed to be at the [-5, 0, 0] location in body fixed coordinates. This is done 
to invoke realism into the algorithm. In a hardware implementation of this algorithm, the 
GPS input would be initially referenced to the antenna location. To perform the lever arm 


correction the INS generated values for the pitch, roll and yaw (©, @®,'¥ ) are used by the 
function "B2NANT.M" to compute the rotation matrix ;R. By multiplying the rotation 


matrix and the position *P,,, and, then, adding the result to the GPS position located at 
the center of gravity of the P-3, results in the following: 


al Rs eee + 3Rx" Payy (351) 
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Figure 3.8: SIMULATOR #1 Block Diagram 


The “SIMULATOR #1” block diagram is shown in Figure 3.8. This block 
diagram performs several functions. First, it performs the integration of INS, GPS and the 
ASM Simulator to derive a target position in the tangent plane coordinate system. Second, 
it performs the same integration again but with true values for the relevant positions and 


Euler angles of the P-3. In this way, the algorithm is validated. Third, this block diagram 


48 


stores the DSTSPI derived target positioning data and the true derived target positioning 


data for post simulation analysis. 
To derive the target position from DSTSPI, the GPS antenna position (* P,,,, ) and 


the Euler angles from INS (©,@®,'¥ ) are used to find the INS/GPS derived position of 
seeker #1 (‘P,, ) in tangent plane coordinates and the INS/GPS derived yaw (‘Y,,,, ) and 
pitch (©,,,,) of the Simulator. This is accomplished by the “P-3 INTEGRATION” 
block as shown in Figure 3.9. Second, this information derived from INS/GPS is then 
used in conjunction with the seecker-to-target information (from the seeker model) to 
derive the target position. The seeker-to-target parameters consist of the range, (Rp,-p ), 
azimuth (AZ — ‘Y,,,,) and elevation (EL —@©,,,, ) between seeker #1 and the target while 


accounting for the attitude of the Simulator. This information is consistent with the 


typical output of a seeker and is discussed in more detail later in this chapter. 

The true target position ("P,,, ), the true P-3 position (* Po, ) and the true Euler 
angles (©',® ,V ), are used to derive the true seeker position (*P,, ), the true Simulator 
yaw (‘Y,,,) and pitch (@,,,) in the “TRUTH MODEL” block. The true position of 
seeker #1, the true yaw and pitch of the Simulator, and the true position of the target are 
used to derive the true seeker-to-target parameters, namely the true range (R ), azimuth 
(AZ —‘¥,,,) and elevation (EL —@,,,,) between seeker #1 and the target. This 


information is used to validate the algorithm in the “TRUE DERIVED TARGET 


POSITION” block. 
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Figure 3.9: P-3 INTEGRATION Block Diagram 


The “P-3 INTERGRATION?” block diagram is shown in Figure 3.9. This block 


performs several functions. First, it takes the GPS position referenced to the GPS antenna 
("P...) and lever arm corrects it to the center of gravity of the P-3 ("Pun ). Next, it 
takes the GPS position referenced to the center of gravity of the P-3 (* P.,.. ) and lever arm 
corrects this position to the forward and after positions of the ASM Simulator #1. These 
positions are denoted as 'P,, and *P,, respectively. After these lever arm corrections 


are made, the forward and after positions of Simulator #1 are known in the tangent plane 


coordinate system. Next, the forward and after positions of Simulator #1 are used to 
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compute the yaw (‘Y,,,, ) and pitch (@,,,, )of the Simulator in the tangent plane coordinate 
system. 
Taking the GPS position at the GPS antenna and lever arm correcting it to the GPS 


position at the center of gravity of the P-3, results in the following: 

"Pee = Pane tpRX°P cen (3-2) 
Where, the value ;Rx*P_..., is computed by the function "B2NCEN.M" using the 
instantaneous values for ©,®,W%. The negative sign indicates that the lever arm 


correction is taken from the antenna to the center of the aircraft. 
Now, by using the GPS position at the center of gravity of the P-3 resolved in 
tangent plane coordinates and lever arm correcting it to the position at seeker #1, the 


following result is obtained: 

SPR = Peet, Rowe, (3-3) 
Where the value }Rx®P., is computed by the function "B2NP2.M" using the 
instantaneous values for ©, ®,'P . 


Likewise, taking the GPS position at the center of gravity of the P-3 and lever arm 


correcting it to the after part of Simulator #1, results in the following: 

/ Pe eR. Rx Po (3-4) 
Where the value }Rx°P,, is computed by the function "B2NP1.M" using the 
instantaneous values for ©,®,%. The functions "B2NP2.M" and "B2NP1.M" are found 


in Appendix A. 
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Given the forward and after positions of the Simulator #1 resolved in tangent plane 
coordinates, the yaw and pitch of Simulator #1 is determined. This is accomplished by 


using the “DETERMINE YAW AND PITCH OF SIMULATOR #1” block which is 


depicted in Figure 3.10. 
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Figure 3.10; DETERMINE YAW AND PITCH OF SIMULATOR #1 Block Diagram 


To determine the yaw and pitch of Simulator #1, first, the difference is taken 
between the forward and after positions of the Simulator. This is given by’ P,,—'P,,- 


This result is a position vector resolved in the tangent plane coordinate system. The length 


of this vector is equivalent to the absolute length of the Simulator. The angles this 


ape 


position vector make with the tangent plane are defined in terms of spherical coordinates 


as shown in Figure 3.11. 


x 
= 





Figure 3.11: Difference between forward and aft Simulator positions resolved in 
spherical coordinates 


From Figure 3.11 we see that the yaw angle is equivalent to the angle 6,,,, and the 
pitch angle is equivalent to@,,,,. The pitch angle is defined to be positive from true 
horizontal to a vertical elevation and negative from true horizontal to straight downward. 
Because the tangent plane coordinate system is defined in a NED orientation, the 
Cartesian to spherical coordinate transformation results in the wrong polarity for the pitch. 
This is why the pitch angle is equal to 0.,,, - 90° vice 90° - 6... The Simulator yaw angle, 


6... 1S defined in terms of the Eulerian notation, ‘P.., Or Voy = 8om - 
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Figure 3.12: TRUTH MODEL Block Diagram 


The block entitled “TRUTH MODEL” is shown in Figure 3.12. This block 
performs the some of the same basic functions as the “P-3 INTEGRATION” block as 


explained earlier in this chapter. First, it takes the true position of the P-3 referenced to 
the center of gravity of the aircraft ("P.,, ) and lever arm corrects this position to the 
forward and after positions of the ASM Simulator #1. These true positions are denoted as 
"P,, and "P,, respectively. Next, it takes forward and aft positions of Simulator #1 and 
computes the yaw (‘Y,,,,) and pitch (@,,,,) of the Simulator in the tangent plane 


coordinates by way of the “GENERATE TRUE AZIMUTH ELEVATION RANGE 
BETWEEN SEEKER AND TARGET®” block as shown in Figure 3.13. Also, the 


difference is taken between the true seeker position and the true target position which 


results in the value 'P,,,—"P,. This value is used in Figure 3.13 to compute the true 
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seeker-to-target parameters namely the range (R_), the azimuth (AZ —‘Y,,,) and the 


elevation (EL — @,,, ). 
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Figure 3.13: GENERATE TRUE AZIMUTH ELEVATION AND RANGE 
BETWEEN SEEKER #1 AND TARGET Block Diagram 
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Now, by using the true position at the center of gravity of the P-3 resolved in tangent 


plane coordinates and lever arm correcting it to the position at seeker #1, the following 


result is obtained: 


>> 


*Pri= PoentaR x Pr, GES) 
Where the value ,;Rx°P,, is computed by the function "B2NP2.M" using the 


instantaneous values for © ,@.. Likewise, taking the true position at the center of 


gravity of the P-3 and lever arm correcting it to the after part of Simulator #1, results in the 


following : 
"Pa= Porn taR x BS (3-6) 
Where the value ,Rx°P,, is computed by the function "B2NP1.M" using the 


instantaneous values for © ,@®.. Next, the true values of the forward and after 
positions of the Simulator #1 and the position vector between seeker #1 and the target are 
used in the SIMULINK block diagram pictured in Figure 3.13 to find the true yaw and 
pitch of the Simulator. To determine the true yaw and pitch of Simulator #1, the 
difference between the forward and after positions of the Simulator, *P,,—'P,,, is used. 
This 1s a position vector resolved in the tangent plane coordinates. Again, the angles this 
position vector make with the tangent plane are defined in terms of spherical coordinates. 
(Refer to Figure 3.11.) Using the same argument that was used in the derivation of the 
INS/GPS derived yaw and pitch, the true yaw angle of Simulator #1 is equivalent to the 
angle .,,, and the true pitch angle of Simulator #1 is equivalent to ©,,,,. These values 
are obtained form the following equations: 


Nae = Ox (3-7) 


Osm = Deni — 90° (3-8) 
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Next, calculate the seeker-to-target parameters by using the value for the difference 


between the true seeker position and the true target position," P,,,—'P,,, and the true yaw 


and pitch of Simulator #1. The value for ™P,,,—'P,, is a position vector resolved in 


tangent plane coordinates. The angles this vector makes with the tangent plane are defined 


by a Cartesian to spherical coordinate transformation. This is seen in Figure 3.14. 





Figure 3.14: Difference between true seeker and target positions resolved in 
spherical coordinates 


Here, the angle AZ is the true azimuth angle between seeker #1 and the target without 
consideration for the yaw angle of the Simulator. The angle EL is the elevation angle 
between seeker #1 and the target without consideration for the pitch angle of the Simulator 


and is defined with respect to horizontal. The values for AZ and EL are given by the 


following equations: 


AZ = Oa (3-9) 


EL = bap — 90° (3-10) 


To find the true azimuth angle between the seeker and the target, the azimuth must 
be compensated by an amount equal to the yaw angle of the Simulator. The goal is to get 
a value for the azimuth angle that is an expected output of the Simulator. A typical radar 
would not account for its orientation with its surroundings when giving its output. The 
result after compensation is AZ —‘¥,,,. The true yaw angle of the Simulator is 
subtracted from the true azimuth angle to maintain consistency with true North. If, for 
example, the Simulator is pointing due East or 90° and the azimuth angle between the 
seeker and the target is also taken to be due East or 90°, then we would expect the 
Simulator to return an angle of zero for the azimuth. 

Now, using the same reasoning, calculate a compensated elevation angle between 
the seeker and the target. This is given as EL —©,,,. 

The true range between the seeker and the target is a easily obtained from the 


Cartesian to spherical coordinate transformation on the position vector "P,,,—'P,, since 


the range does not care about what the Simulator is doing. This is givenas R . 

To validate the algorithm, the “TRUE DERIVED TARGET POSITION” block 
diagram as pictured in Figure 3.15 is used. The goal of this block diagram is to find the 
target position in tangent plane coordinates using true values for inputs. This is 


accomplished be first taking the true seeker #1 to target parameters and correcting these 
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values for the true yaw and pitch of the Simulator. This allows for the extraction of the 


true position vector between the seeker and the target via a coordinate transformation. 
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Figure 3.15: TRUE DERIVED TARGET POSITION Block Diagram 


Given the true seeker to target parameters, R ,AZ —‘¥,,, and EL —©,,,, first 
these values must be corrected for the true yaw and pitch of the Simulator. This is done by 
adding the respective values in the following manner: 

AZ =AZ -¥,,, + Yom (3-11) 
EL = EL -0,,,+9om (3-12) 
If this correction were not made, the spherical to Cartesian coordinate transformation 


would result in severe displacement errors from the yaw and pitch of the Simulator. The 


elevation angle, as mentioned earlier, is defined with respect to honzontal. To put the 


elevation term in spherical coordinates simply add 90 degrees, or $,,, = 90° +EL. 
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Next, perform a spherical to Cartesian coordinate transformation. The result in 
tangent plane coordinates is the position vector "P,,.—'P,,. By adding the true position 
of seeker #1 to this quantity, the result is the true target position in tangent plane 
coordinates. This is seen in the following equation: 

[Pe P= bee Pe (3-13) 

The next step in the algorithm development is to talk about the seeker model for 
the conical scan radar. The block diagram entitled “SEEKER MODEL FOR CONICAL 
SCAN RADAR?” is shown in Figure 3.16. The three parameters that are relevant to this 
algorithm are the range, elevation and azimuth angle between the modeled seeker and the 
target while accounting for the ASM Simulator's attitude. The output of a real seeker, as 
mentioned earlier, does not account for its surroundings. It only knows the target in a 


relative sense. 
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Figure 3.16: SEEKER MODEL FOR CONICAL SCAN RADAR Block Diagram 


Because the information received from the seeker is not referenced to a particular 
coordinate system, true values for the position of the seeker are used to derive the seeker 


to target parameters. Also, because the attitude of the seeker is not accounted for in its 


output, true values for the yaw and pitch of the seeker are used. 


Referring to Figure 3.16, first the difference is taken between the true seeker 


position and the true target position which results in the value "P,,,—'P,,. This value is 


used to help compute the modeled seeker to target parameters. The value for *P,,,—'P,,, 


is a position vector resolved in tangent plane coordinates. By performing a Cartesian to 
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spherical coordinate transformation on this position vector, the values R, $,,,, and 


O.,4z are obtained. Note, that these values are based on the true positions of the seeker 
and the target as denoted by the prime. 


Next, the value for Ris quantized in 100 meter segments to simulate the range 
resolution of the modeled seeker. The output of the quantizer is R,,,,, which is the range 
between the seeker and the target that has an error component based on the range 
resolution of the seeker. This quantity is defined as 

Roser = Reare (3-14) 
where “DER” is defined as a derived value. 

The MATLAB function "CONTHETA.M" adjusts the precision in the angular 
resolution of the seeker based on the signal to noise ratio and the scintillation as discussed 
in Chapter II. A similar MATLAB function, “THETAD.M,” performs the same function 
as ““CONTHEATA.M” but for the monopulse seeker model found in the block diagram 
entitled “SIMULATOR #2” The signal to noise ratio changes as a function of the range 
between the seeker and the target. The angular resolution is normally distributed resulting 
in the value 60. The value for 60 contributes errors in both directions of the solid angle 
subtended by the radar beam and, therefore, must be added to both 0,,, and Or,p - 

The angle AZ is the azimuth angle between seeker #1 and the target without 
considering the yaw angle of the Simulator. Likewise, the angle EL is the elevation angle 
between seeker #1 and the target without considering the pitch angle of the Simulator and 
is defined with respect to horizontal. The values for AZ and EL are given by the 


following equations: 
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AZ = 0.,,, +60 (3-15) 
EL = 67,2 +50 — 90° (3-16) 

The azimuth angle between the seeker and the target after compensating for the 
true yaw of the Simulator is AZ—‘Y,,,,. Again, the true yaw angle of the Simulator is 
subtracted from the azimuth angle to maintain consistency with true North. The elevation 
angle between the seeker and the target after compensating for the pitch angle of the 
seeker is given by EL-©,,,,. 

The final step in the algorithm is derive the target position from the INS/GPS 
derived parameters and the seeker to target parameters. Shown in Figure 3.17 is the block 
diagram entitled “INS/GPS/RADAR DERIVED TARGET POSITION”. This block 
performs the same basic functions as Figure 3.15 with a few minor modifications. 

Given the derived seeker-to-target parameters, Rppp, AZ — Vy, and EL-®,,,; 
these values must be corrected for the INS/GPS derived yaw and pitch of the Simulator. 
This is done by adding the respective values in the following manner: 

AZ pep = AZ — Yom + Yon = Oper (3-17) 

EL per = EL—- Og, + Osim (3-18) 
Notice that yaw and pitch terms of the Simulator do not fully cancel. The residual values 
are due to the errors involved with INS, GPS and the seeker. These errors translate 
directly into displacement errors in fixing the target’s position in tangent plane 


coordinates. These errors are discussed in more detail in the following chapter entitled 


“VALIDATION AND ANALYSIS.” 
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Figure 3.17: INS/GPS/RADAR #1 DERIVED TARGET POSITION Block Diagram 


The elevation term, ELp,,, 1S not in the proper form. To put the elevation term in 
the form of spherical coordinates simply add 90 degrees, or 9p.2 = 90°+EL,,,. This 
results in all values in terms of spherical coordinates. 

Next, perform a spherical to Cartesian coordinate transformation. The result in 
tangent plane coordinates is the position vector "P,,..—'P,,. To derive the target position 
in tangent plane coordinates, add the INS/GPS derived position of seeker #1 to this 


quantity. This is seen in the following equation: 


*Proer= Pracr— Pr +P, G-y 
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Note, the true and derived position vectors of seeker #1 do not cancel. This is due 
to the fact that the derived seeker-to-target parameters are based on the true seeker 
position. 

Next, assume the target has an upper limit on its speed. This means the x, y and z 
components of the DSTSPI derived target position can be rate limited to a specified value. 
In other words, the target position is not allowed to change by more than a specified 
amount. This rate of change is limited to a value of 50 KTS. The z component of the 
target position could have been limited to values close to the surface, but, this would have 


possibly resulted in the miss detection of any airborne decoys such as a chaff cloud. 
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IV. VALIDATION AND ANALYSIS 


All the figures in this chapter were generated using MATLAB. Also, all the 


figures were generated under the following conditions: 


The P-3’s initial course is 000° true 

The P-3’s initial altitude is 2000 feet 

The P-3’s initial speed is 500 ft/sec or 296.242 kts 

The P-3’s initial position in the tangent plane coordinate system is [0, 0,-2000] 
Mpicer 

The Target’s initial course in Chapter IV, Section’s A through C is 270° true 
The Target’s speed is 22 kts throughout the simulation 

The Target’s initial location in the tangent plane coordinate system is [5, 0, 0] 
in nautical miles 


All data collection periods were 60 seconds 


A. INS AND GPS ERRORS 


This section looks at the modeled errors from both INS and GPS. Plots of these 


errors were generated using MATLAB. As discussed in Chapter II, there are several 


factors which contribute to the errors seen from INS and GPS. One of the limiting factors 


in the algorithm is the fact that these errors are always present. The biggest contributors to 


the displacement error in the derived target position are due to the errors from INS and the 


seekers. Although GPS errors are important, an error in the seeker's position only results 


in the same order magnitude error in the derived target's position. Simply stated, because 


the arc length is equal to the angle in radians multiplied by the radius (or the distance 
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between target and the seeker), any angular errors prove to be detrimental to the algorithm. 
These angular errors result from both the seeker’s ability to measure the azimuth and 
elevation between the Simulator and the target and the INS’s ability to measure the Euler 
angles. As a hypothetical case, assume that a particular seeker's position is known to 
within one meter and the pitch, roll and yaw of the Simulator are known to within 0.5°, or 
0.009 radians. If the seeker were 5 NMI from the target, this would translate to a cross 
range displacement error in the derived target position of roughly 83 meters assuming the 
seeker measured the azimuth and elevation precisely. 

As mentioned in Chapter II, DGPS, theoretically, should provide positions to 
within a one meter accuracy. Figures 4.1 through 4.3 show the GPS position errors 
relative to a NED orientation. From these Figures one can conclude the maximum 
amplitude for the displacement error occurs at roughly 12 seconds into the simulation over 
the sample period. At this point the value for this amplitude is about 2.5 meters. These 


errors are roughly consistent with actual DGPS performance specifications. 
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Figure 4.1: X-component of the GPS displacement error 
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Figure 4.2: Y-component of the GPS displacement error 
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Figure 4.3: Z-component of the GPS displacement error 


The errors in the INS values for pitch roll and yaw are depicted in Figures 4.4 
through 4.6. From Figures 4.4 through 4.6, the maximum error in the INS’s computation 
of the Euler angles is approximately +/- 0.4 degrees. This equates to a cross range 
displacement error in the derived target position of approximately 65 meters at a 5 NMI 
range. The errors in the Euler angle computation for this model are comparable to the 
errors found from an actual INS system. Specifically, NRL’s P-3 uses the LTN-72 which 
is an INS with an advertised maximum amplitude for the Euler angle measurement error of 


less than 0.5 degrees. [Ref. 8, p. 35] 
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Figure 4.4: INS pitch error of the P-3 
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Figure 4.5: INS roll error of the P-3 
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Figure 4.6: INS yaw error of the P-3 


B. SEEKER POSITION ERRORS 


Figures 4.7 through 4.12 show the position errors of seeker #1 and seeker #2. 
Notice, that these curves follow the errors associated with GPS very closely with slight 
variations due to the INS errors. As mentioned earlier, these errors translate to the same 
order magnitude errors in the derived target positions. There are two contributions to the 
overall error in deriving the position of each seeker. The first is the contribution from 
DGPS. The uncertainty in the DGPS position at the GPS antenna is lever arm corrected to 
the seeker position. In other words, the position errors of DGPS translate directly into 


position errors of the seekers. 


gos 


Additionally, any angular errors from INS contribute to the overall position error 
of each seeker. A 0.5 degree error in any one of the Euler angle computations results in a 
5 centimeter error in the position of each seeker for a lever arm of 5 meters. This means 
that the errors from INS are insignificant in deriving each seeker’s position. This does not 
mean that the Euler angles are themselves are insignificant. Without the lever arm 
correction, the errors in deriving the seeker’s position in a fixed inertial coordinate system 


would be substantially greater. 
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Figure 4.7: Seeker #1 position error in North/South direction 
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Figure 4.8: Seeker #1 position error in East/West direction 
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Figure 4.9: Seeker #1 position error in Up/Down direction 
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Figure 4.10: Seeker #2 position error in North/South direction 
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Figure 4.11: Seeker #2 position error in East/West direction 
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Figure 4.12: Seeker #2 position error in Up/Down direction 


Figures 4.13 and 4.14 provide a comparison of the true and derived positions of 
seeker #1. In Figures 4.13 and 4.14 the initial displacement in the +y direction for the 
position of seeker #1 is due to the location of seeker #1 in body fixed coordinates relative 
to the center of the P-3. Additionally, the P-3 does not fly a straight line. The course 
variations are due to the “P-3 FLIGHT PROFILER” model. Because this model 
continuously tries to reduce the error between the ordered course and the actual course, the 
P-3 can not maintain a steady course but continuously hunts for the ordered course. 

Figures 4.15 and 4.16 provide a comparison of the true and derived positions of 
seeker #2. The initial displacement in the -y direction is due to the placement of seeker #2 


with in the body fixed coordinates. 
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Figure 4.13: True position of seeker #1 
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Figure 4.14: INS/GPS derived position of seeker #1 
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Figure 4.15: True position of seeker #2 
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Figure 4.16: INS/GPS derived position of seeker #2 
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C. SEEKER CONTRIBUTED ERRORS 


Figure 4.17 shows the error from the range resolution for the conical scan seeker 
model used in the simulation. Because the range resolution is modeled as 100 meters, an 
envelope of 50 meters is evident throughout the simulation. Notice that the range error 
does not decrease over the sample interval. This 1s due to the fact that the range resolution 
is a function of the seeker hardware and not the signal to noise ratio. The range resolution 
of this seeker model translates to down range errors 1n the derived target position as 1s seen 


in the following section. 
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Figure 4.17: Range error from conical scan seeker (COSRO) 


Figure 4.18 shows the precision in the measurement of azimuth and elevation 
angles from the conical scan radar. This precision is approximated as +/- 0.15 degrees 


over this data collection period. This precision in the measurement of the azimuth and 
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elevation remains relatively constant during this period because the aircraft is closing the 
target from 5 NMI. Because the signal to noise ratio diminishes the effects of the error 
slope and crossover loss in the conical scan seeker, at this range, the errors in the azimuth 
and elevation are dominated by the scintillation error. 

An angular precision of +/- 0.15 degrees is expected to result in cross-range errors 


on the order of 25 meters when the seeker to target separation is 5 NMI. 


AZIMUTH AND ELEVATION ERROR (degrees) 
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Figure 4.18: Azimuth and elevation error from conical scan seeker (COSRO) 


Figure 4.19 shows the error in the range measurement for the monopulse seeker or 
seeker #2. Here the range resolution is 20 meters. As expected, the magnitude of the 
range error is +/- 10 meters. Notice the crossover rate for the range error is 5 times faster 
than for the crossover rate of the conical scan range error. This results from having a 


range resolution that is 5 times smaller in the case of the monopulse seeker. Again, the 
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range resolution of the monopulse seeker results in down range errors in the derived target 


position of roughly the same order magnitude. 
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Figure 4.19: Range error from monopulse seeker 


Figure 4.20 shows the azimuth and elevation measurement errors for the 
monopulse seeker. Because the monopulse seeker is not affected by scintillation errors, an 
increase in signal to noise ratio directly results in an increase in the precision of the 
measurement. As the range between the seeker and target decreases, the azimuth and 
elevation measurement errors decrease from +/- 0.05 degrees to a negligible amount. At 
0.05 degrees the cross range error in the derived target position would be roughly 8 meters 


with a seeker to target separation of 5 NMI. 
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Figure 4.20: Azimuth and elevation error from monopulse seeker 


D. DERIVED TARGET POSITION ERRORS 


This section shows the true and derived target positions for each seeker type in the 
tangent plane coordinate system for several different types of target motion. In addition, 
this section shows the position errors that result from the DSTSPI algorithms for each 
case. 

Figures 4.21 through 4.29 were generated with the target dead in the water. 
Figures 4.30 through 4.38 were generated with the target on a course of 270° true. Figures 
4.39 through 4.47 were generated with the target on a course of 090° true. Figures 4.48 
through 4.56 were generated with the target on a course of true North. Figures 4.57 
through 4.65 were generated with the target on a course of 180° true. Figures 4.66 through 


4.74 were generated with the target making a slow turn to port from 000° to 270° true. 
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Figures 4.75 through 4.83 were generated with the target making a slow turn to starboard 
from 000° to 090° true. Figures 4.84 through 4.93 were generated with the target making 
a slow turn to starboard from 180° to 270° true. Figures 4.94 through 4.101 were 
generated with the target making a slow turn to port from 180° to 090° true. 


Table 4.1 is provided as a summary of the errors from the DSTSPI algorithms for 


seeker #1 and seeker #2 for each of the nine types of target motion. 
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Table 4.1: Summary of DSTSPI derived target position errors 
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From Table 4.1 several conclusions are made. The errors between the derived 


target positions from both seekers and the true target position are essentially constant over 
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a wide range of target motions. The errors are greatest in the cross range direction. The 
errors from the COSRO seeker are larger than for the monopulse seeker. This is especially 
true when the seeker-to-target separation is small. 

Assume that the errors that from the DSTSPI algorithms are roughly the same as 
the maximum values seen in Table 4.1. If this assumption is used, then the area of 
uncertainty in deriving a target’s position is roughly a sphere of radius 40 meters around 
the true target position. If the target’s speed is 22 kts or 11.3178 meters/sec, then the 
elapsed time before any decoys are detected after launch is about 3.5 seconds. This 
latency in detecting the decoy is defined as the time resolution, dt, of the DSTSPI 
algorithms. The time resolution, dt, decreases as the area of uncertainty of the target’s 
position decreases. From the figures presented in this section, the general trend of the area 
of uncertainty of the target’s position is that it decreases with a decreasing seeker-to-target 


separation. 
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Figure 4.21: True target position - DIW 
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Figure 4.22: Simulator #1 DSTSPI derived target position - DIW 
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Figure 4.23: Simulator #2 DSTSPI derived target position - DIW 
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Figure 4.24: Simulator #1derived target position error in the North/South direction - 
DIW 
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Figure 4.25: Simulator #1 derived target position error in the East/West direction - 
DIW 


Z COMPONENT OF ERROR (meters) 





10 20 30 40 50 60 


Time (second) 


Figure 4.26: Simulator #1 derived target position error in the Up/Down direction - 
DIW 
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Figure 4.27: Simulator #2 derived target position error in the North/South direction 
- DIW 
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Figure 4.28: Simulator #2 derived target position error in the East/West direction - 
DIW 
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Figure 4.29: Simulator #2 derived target position error in the Up/Down direction - 
DIW 
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Figure 4.30: True target position - CSE 270° 
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Figure 4.31: Simulator #1 DSTSPI derived target position - CSE 270° 
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Figure 4.32: Simulator #2 DSTSPI derived target position - CSE 270° 
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Figure 4.33: Simulator #lderived target position error in the North/South direction - 
CSE 270° 
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Figure 4.34: Simulator #1 derived target position error in the East/West direction - 
CSE 270° 
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Figure 4.35: Simulator #1 derived target position error in the Up/Down direction - 


CSE 270° 
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Figure 4.36: Simulator #2 derived target position error in the North/South direction 
- CSE 270° 
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Figure 4.39: True target position - CSE 090° 
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Figure 4.40: Simulator #1 DSTSPI derived target position - CSE 090° 
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Figure 4.41: Simulator #2 DSTSPI derived target position - CSE 090° 
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Figure 4.42: Simulator #1 derived target position error in the North/South direction - 
CSE 090° 
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Figure 4.43: Simulator #1 derived target position error in the East/West direction - 
CSE 090° 
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Figure 4.44: Simulator #1 derived target position error in the Up/Down direction - 
CSE 090° 
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Figure 4.45: Simulator #2 derived target position error in the North/South direction 
- CSE 090° 
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Figure 4.46: Simulator #2 derived target position error in the East/West direction - 
CSE 090° 
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Figure 4.48: True target position - CSE 000° 
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Figure 4.49: Simulator #1 DSTSPI derived target position - CSE 000° 
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Figure 4.50: Simulator #2 DSTSPI derived target position - CSE 000° 
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Figure 4.51: Simulator #1lderived target position error in the North/South direction - 
CSE 000° 
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Figure 4.52: Simulator #1 derived target position error in the East/West direction - 
CSE 000° 
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Figure 4.53: Simulator #1 derived target position error in the Up/Down direction - 
CSE 000° 
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Figure 4.54: Simulator #2 derived target position error in the North/South direction 
- CSE 000° 
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Figure 4.55: Simulator #2 derived target position error in the East/West direction- 
CSE 000° 
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Figure 4.56: Simulator #2 derived target position error in the Up/Down direction - 
CSE 000° 
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Figure 4.57: True target position - CSE 180° 
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Figure 4.58: Simulator #1 DSTSPI derived target position - CSE 180° 
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Figure 4.59: Simulator #2 DSTSPI derived target position - CSE 180° 
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Figure 4.60: Simulator #1derived target position error in the North/South direction - 
CSE 180° 
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Figure 4.61: Simulator #1 derived target position error in the East/West direction - 
CSE 180° 
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Figure 4.62: Simulator #1 derived target position error in the Up/Down direction - 
CSE 180° 
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Figure 4.63: Simulator #2 derived target position error in the North/South direction 
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Figure 4.65: Simulator #2 derived target position error in the Up/Down direction - 
CSE 180° 
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Figure 4.66: True target position - CSE 000° To 270° 
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Figure 4.67: Simulator #1 DSTSPI derived target position - CSE 000° To 270° 
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Figure 4.68: Simulator #2 DSTSPI derived target position - CSE 000° To 270° 
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Figure 4.69: Simulator #1 derived target position error in the North/South direction - 
CSE 000° To 270° 
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Figure 4.70: Simulator #1 derived target position error in the East/West direction - 
CSE 000° To 270° 
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Figure 4.71: Simulator #1 derived target position error in the Up/Down direction - 
CSE 000° To 270° 
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Figure 4.72: Simulator #2 derived target position error in the North/South direction 
- CSE 000° To 270° 
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Figure 4.73: Simulator #2 derived target position error in the East/West direction - 
CSE 000° To 270° 
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Figure 4.74: Simulator #2 derived target position error in the Up/Down direction - 
CSE 000° To 270° 
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Figure 4.75: True target position - CSE 000° To 090° 
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Figure 4.76: Simulator #1 DSTSPI derived target position - CSE 000° To 090° 
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Figure 4.77: Simulator #2 DSTSPI derived target position CSE 000° To - 090° 
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Figure 4.78: Simulator #1derived target position error in the North/South direction - 
CSE 000° To 090° 
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Figure 4.79: Simulator #1 derived target position error in the East/West direction - 
CSE 000° To 090° 
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Figure 4.80: Simulator #1 derived target position error in the Up/Down direction - 
CSE 000° To 090° 
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Figure 4.81: Simulator #2 derived target position error in the North/South direction 
- CSE 000° To 090° 
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Figure 4.82: Simulator #2 derived target position error in the East/West direction - 
CSE 000° To 090° 
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Simulator #2 derived target position error in the Up/Down direction - 
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Figure 4.84: True target position - CSE 180° To 270° 
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Figure 4.85: Simulator #1 DSTSPI derived target position - CSE 180° To 270° 
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Figure 4.86: Simulator #2 DSTSPI derived target position - CSE 180° To 270° 
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Figure 4.87: Simulator #1derived target position error in the North/South direction - 
CSE 180° To 270° 
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Figure 4.88: Simulator #1 derived target position error in the East/West direction - 
CSE 180° To 270° 
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Figure 4.89: Simulator #1 derived target position error in the Up/Down direction - 
CSE 180° To 270° 
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Figure 4.90: Simulator #2 derived target position error in the North/South direction 
- CSE 180° To 270° 
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Figure 4.91: Simulator #2 derived target position error in the East/West direction - 
CSE 180° To 270° 
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Figure 4.92: Simulator #2 derived target position error in the Up/Down direction - 
CSE 180° To 270° 
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Figure 4.93: True target position - CSE 180° To 090° 
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Figure 4.94: Simulator #1 DSTSPI derived target position - CSE 180° To 090° 
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Figure 4.95: Simulator #2 DSTSPI derived target position - CSE 180° To 090° 
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Figure 4.96: Simulator #1derived target position error in the North/South direction - 
CSE 180° To 090° 
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Figure 4.97: Simulator #1 derived target position error in the East/West direction - 
CSE 180° To 090° 
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Figure 4.98: Simulator #1 derived target position error in the Up/Down direction - 
CSE 180° To 090° 
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Figure 4.99: Simulator #2 derived target position error in the North/South direction 
- CSE 180° To 090° 
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Figure 4.100: Simulator #2 derived target position error in the East/West direction - 
CSE 180° To 090° 
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Figure 4.101: Simulator #2 derived target position error in the Up/Down direction - 
CSE 180° To 090° 
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V. HARDWARE IMPLEMENTATION 


As discussed in Chapter III, the DSTSPI algorithms require specific inputs. These 
inputs are a GPS position relative to the location of the GPS antenna, the pitch, roll and 
yaw of the aircraft that 1s given by the INS, and the azimuth, elevation and range between 
the seeker and the target. These inputs are shown in Figure 5.1. Of the three inputs, only 


the GPS position presents any real problems. 
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Figure 5.1: Hardware implementation of the DSTSPI algorithms 
The algorithms presented in Chapter III assumed all fixed inertial positions are 
given in tangent plane coordinates. This is done as a matter of convenience. To continue 
to use the tangent plane coordinates in the hardware implementation, a fixed position 
located in geodetic coordinates must be defined somewhere in the ASW/EW engagement 
test scenario. This fixed point of reference could be located at the center of the operating 


area that is allocated prior to the test. This point is taken from the existing navigational 


WA 


chart, the coordinates of which are given in the geodetic sense of Latitude and Longitude. 
The actual location is arbitrary. 

To locate the P-3 in tangent plane coordinates we must take the position of the 
origin of the tangent plane coordinate system given in geodetic coordinates and perform a 
coordinate transformation from geodetic to ECEF coordinates. This position now 
becomes the origin of the tangent plane coordinate system given in ECEF coordinates. 

The geodetic to ECEF coordinate transformation is easily computed. First, three 
quantities must be defined. The first of these quantities is the flattening factor, f, which 
represents the relative flatness of the reference ellipsoid. The mathematical definition for f 


is defined as (Ref. 1, p. 118-119] 





f= (5-1) 


where a and b are the semi-major and semi-minor axes of the reference ellipsoid, 
respectively. The next quantity that must be defined is the eccentricity, ¢. The 
eccentricity of the reference ellipsoid is defined as 

se =2f—f’. (5-2) 
The last quantity that must be defined is the normal, N, which represents the length of the 


ellipsoidal normal from the ellipsoidal surface to the intersection of the polar axis (Zgcpr). 


Mathematically, N is defined as 


a 


— 5-3 
yl—e’ sin’ > (ag) 


where @ is the geodetic latitude. Now, by using these three quantities, the geodetic to 


N= 


ECEF transformation is defined by the following equations 
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Xecer = (N + h) cos cosa (5-4) 

Yecer =(N +h) cos sina (5-5) 

Zecer =|N(1-e7)+h|sing (5-6) 
where h is the geodetic height or height above reference ellipsoid and A is the geodetic 
longitude. Equations (5-4) through (5-6) are incorporated into the MATLAB function 
“LL2ECEF.M” found in Appendix A. 

Now, the next step is to find the position of the P-3 relative to the tangent plane 
coordinates. To do this we take the GPS position given in ECEF coordinate and subtract it 
from the origin of the tangent plane given in ECEF coordinates. This results in a position 
vector that has the correct orientation, NED. 

Another option is to use a DGPS reference station. Here a portable DGPS 
reference station 1s located as near as possible to the testing location (this location must be 
on land). This station provides the necessary differential corrections to the local receiver. 
The advantage of using a portable DGPS reference station is that they usually provide 
better corrections than the Coast Guard differential stations. Another advantage is that 
these remote stations have the ability to record the GPS data. This would allow for further 
reductions in the errors from GPS. [Ref. 3] A typical system is found in Appendix B. The 
problem is that these reference stations must be setup and in operation 24 hours prior to 
the test scenario. Another problem is the cost. A typical reference station costs 
approximately twenty thousand dollars. 

Since the algorithms presented in Chapter III are not intended for real time use, an 


actual hardware implementation only involves recording the three inputs for post-test 
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analysis. The difficulty in doing this is to ensure the data is saved in a format that is 
compatible with a particular post analysis software. To make the process easier, several 
GPS receiver manufacturers provide software packages that provide a post analysis 
capability. To this end it behooves the user to pick a GPS system that comes with a 
complete software package. 

Another potential hazard in the recording of the data for post-test analysis is to 
ensure that all the recorded times are properly synchronized. Any clock differences 
translate directly into position errors. This is especially true for the time-tagged data of 
the INS and the Simulators since the angular rates can be quite large. 

Perhaps the most difficult problem is the determination of the actual target position 
as a function of time given noisy data from the DSTSPI algorithms. The goal of the post- 
test analysis 1s to provide sufficient accuracy in the determination of the target’s position 
so as to allow for the recognition of any launched decoys. To help attain this goal it is 
possible to take the raw data from the DSTSPI algorithms and in some way make an 
educated guess about the true target position. A first order approximation in obtaining the 
target’s true position might be a simple least-squares fit of the data, but because of the 
cross range errors which are a consequence of the INS and seeker errors, a simple least- 
squares fit on the data does not provide any useful information. This is especially true for 
a maneuvering target. A better method in determining the target’s true position given the 
noisy data might involve an iterative process by which the errors in all directions are 


minimized at the same rate. 
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VI. CONCLUDING REMARKS 


The effectiveness of the at sea Anti-Ship Missile/Electronic Warfare engagement 
test scenario is heavily dependent upon the results of this thesis. By knowing the target's 
position in a fixed inertial coordinate system, the post analysis problem of distinguishing 
between target and decoy is made easier. Any launched decoys during the test scenario 
show up as an instantaneous deviation in the target track. This allows the test engineers 
determine if the seeker broke lock on the target and tracked an offboard decoy (e.g., chaff). 
By improving the method by which the at sea ASM/EW test scenario evaluation is 
performed, the war-fighting capabilities of the U. S. Navy improve. 

Also, by knowing a target's position relative to a fixed inertial coordinate system, a 
target's position can be sent to any platform and be understood immediately without being 
specific to the platform’s Fire Control System. Target positions that are fixed in a fixed 
inertial coordinate system allow for ease of processing because the information is sensor 
independent. 

Considering the results obtained in this simulation, it is feasible to derive a target’s 
position in a fixed inertial coordinate system from the integration of INS, GPS and radar. 
The precision to which the target’s position is derived depends heavily upon the quality of 
the navigation equipment used in the integration process. 

With the advent of more and more sophisticated weaponry, it becomes increasingly 
more important to provide the best available means by which an evaluation of these 
weapons is made. Of particular significance is the foreign missile threat. It 1s here that a 


great deal of time and money is spent in determining the effectiveness of our shipboard 
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self defense systems against these threats. The validation of a shipboard self defense 


system against a missile threat usually involves modeling and simulation. 
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APPENDIX A 


MATLAB FILES: 


AZIMUTH.M 


function out = azimuth(theta): 

%ANDREW ROWE 

%THIS FUNCTION COMPUTES THE AZIMUTH ANGLE 
%RELATIVE TO TRUE NORTH FROM THE THETA GIVEN 
%FROM CARTESIAN TO SPHERICAL COORDINATE 
%TRANSFORMATION 

%10/10/96 

x = theta; 

y=450-x; 

if y >= 360 

bay 000; 

end; 

out = y; 


B2NANT.M 


function out = b2nant(vec): 

% ANDREW ROWE 

% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 

% ROLL=PHI PITCH=THETA YAW=PSI 

% pl IS THE POSITION VECTOR OF THE GPS ANTENNA IN BODY FIXED FRAME 
% COORDINATES RELATIVE TO THE CENTER OF THE AIRCRAFT 

% 10/10/96 

phi = vec(1); 

theta = vec(2); 

psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)]; 

Ttheta = [cos(theta) 0 sin(theta); 0 1 0; -sin(theta) 0 cos(theta)]; 

Tpsi= [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1); 

p1=[-5;0;0]; 

out = inv(Tpsi*Ttheta*Tphi)*p1; 


B2NCEN.M 


function out = b2ncen(vec): 

% ANDREW ROWE 

% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 

% pl IS THE POSITION VECTOR OF THE CENTER OF THE AIRCRAFT RELATIVE TO THE 

% POSITION OF THE GPS ANTENNA 

% ROLL=PHI PITCH=THETA YAW=PSI 

% 10/10/96 

phi = vec(1); 

theta = vec(2); 
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psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)]; 

Ttheta = (cos(theta) 0 sin(theta); 0 1 0; -sin(theta) 0 cos(theta)]; 
Tpsi = [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1); 
p1=[5;0;0]; 

out = inv(Tpsi*Ttheta*Tphi)*p1; 


B2ZNP1.M 


function out = b2np1(vec); 

% ANDREW ROWE 

% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 

% pl IS THE POSITION VECTOR OF THE AFT PART OF THE SIMULATOR #1 RELATIVE TO 
% THE CENTER OF THE AIRCRAFT IN BODY COORDINATES 

% ROLL=PHI PITCH=THETA YAW=PSI 

% 10/10/96 

phi = vec(1); 

theta = vec(2); 

psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)]; 

Ttheta = [cos(theta) 0 sin(theta); 0 1 0; -sin(theta) 0 cos(theta)]; 

Tpsi = [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1); 

p1=[-2;5;0]; 

out = inv(Tpsi*Ttheta*Tphi)*p1; 


B2NP12.M 


function out = b2np12(vec); 

% ANDREW ROWE 

% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 

% pl IS THE POSITION VECTOR OF THE AFT PART OF THE SIMULATOR #2 RELATIVE TO 
% THE CENTER OF THE AIRCRAFT IN BODY COORDINATES 

% ROLL=PHI PITCH=THETA YAW=PSI 

% 10/10/96 

phi = vec(1); 

theta = vec(2); 

psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)]; 

Ttheta = [cos(theta) 0 sin(theta); 0 1 0; -sin(theta) 0 cos(theta)]; 

Tpsi = [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1); 

p1=[-2;-5;0]; 

out = inv(Tpsi*Ttheta* Tphi)*p1; 


B2NP2.M 


function out = b2np2(vec): 
% ANDREW ROWE 


% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 

% pl IS THE POSITION VECTOR OF THE FORWARD PART OF THE SIMULATOR #1 

% RELATIVE TO THE CENTER OF THE AIRCRAFT IN BODY COORDINATES 
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% ROLL=PHI PITCH=THETA YAW=PSI 

% 10/10/96 

phi = vec(1); 

theta = vec(2); 

psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)}]; 

Ttheta = [cos(theta) 0 sin(theta); 0 1 0; -sin(theta) 0 cos(theta)]; 
Tpsi = [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1); 
ml=|2;5;0); 

out = inv(Tpsi*Ttheta*Tphi)*p1; 


B2NP22.M 


function out = b2np22(vec): 

% ANDREW ROWE 

% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 

% pl IS THE POSITION VECTOR OF THE FORWARD PART OF THE SIMULATOR #2 
% RELATIVE TO THE CENTER OF THE AIRCRAFT IN BODY COORDINATES 

% ROLL=PHI PITCH=THETA YAW=PSI 

% 10/10/96 

phi = vec(1); 

theta = vec(2); 

psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)]; 

Ttheta = [cos(theta) 0 sin(theta); 0 1 0; -sin(theta) 0 cos(theta)]; 

Tpsi = [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1]; 

p1=—(2;-5;0); 

out = inv(Tpsi* Ttheta*Tphi)*p1; 


CONTHETA.M 


function out = contheta(component); 

% ANDREW ROWE 

% THIS FUNCTION COMPUTES THE ANGULAR PRECISION FOR A CONICAL SCAN RADAR 
% SEEKER. THIS IS USED IN THE SEEKER MODEL 

% 10/10/96 

R = component(1); “range in meters 


rae 250e3: %power in watts 250 Kw 
sigma = 3000; Yoradar cross section square meters 
f= 9e9; Yradar frequency of 9 Ghz 
alpha = .055; “attenuation dB/km page279 Barton 
Gain = 1995.26; % antenna gain 33 dB 
Cc = 2.99798; Y%speed of light 
lambda = c/f; Yowavelenght in meters 
B = 10e6; Ynoise bandwidth 
P— 10; %noise figure 
kT = 4.14e-21; Yjoules at 300 K 
bw = .0785398; %3dB transmitter beam width (uniform illumination) 4.5 degrees 
SN = Pt*Gain* Gain* sigma* 10“(-0.2*alpha* R/1000)*lambda*lambda/4/pi/kT/B/F/4/4/pi/p/R/R/R/R; 
sigma s=0.0125*bw; “scintillation error at optimized squint angle 
“Barton page 390 
SN = 100*SN; % 100 pulse integration improvement 


ey 


thetasig = bw/(SN)*0.5 + sigma _s; Yoradians 
THETAD.M 


function thetasig = thetad(component); 

% ANDREW ROWE 

% THIS FUNCTION COMPUTES THE ANGULAR PRECISION FOR A MONOPULSE 
% SEEKER. THIS IS USED IN THE SEEKER MODEL 

% 10/10/96 


R = component(1); %range in meters 

Pt = 30e3; “power in watts 30 Kw 

sigma = 3000; “radar cross section square meters 

f= 17e9; Yradar frequency of 17 Ghz 

alpha = .055; “attenuation dB/km page279 Barton 

Gain = 218.77; “% antenna gain 23.4 dB 

C=Z29979e8: “speed of light 

lambda = c/f; “wavelenght in meters 

B = 10e6; %noise bandwidth 

|e rele “noise figure 

kT = 4.14e-21; %joules at 300 K 

bw = .13962; %3dB transmitter beam width (uniform illumination) 8 degrees 
SN = Pt*Gain*Gain*sigma* 10“(-0.2*alpha*R/1000)*lambda*lambda/4/pi/kT/B/F/4/4/pi/p/R/R/RIR; 
km = 1.9; Yetror slope for monopulse radar dual horn 

thetasig = bw/km/(100*2*SN)‘0.5; Y%radians 


HEAD.M 


function out = head(direc): 
% ANDREW ROWE 


% THIS FUNCTION DETERMINES THE HEADING BASED ON NOT EXCEEDING 180 DEGREE 
% TURN DURING AIRCRAFT MANEUVERS 
% 10/10/96 
x = direc(1); 
y = direc(2); 
if (x > y) & ((x-y) < pi), 
direction = 1; 
elseif (x < y) & ((y-x) < pi), 
direction = -1; 
elseif (x > y) & ((x-y) >= pi), 
direction = -1; 
elseif (x < y) & ((y-x) >= pi), 
direction = 1; 
elseif x==y, 
direction = 1; 
end; 
out = direction; 


OVER2PI.M 


function out = over2pi(direc): 
% ANDREW ROWE 


% THIS FUNCTION ENSURES THAT THE HEADING FALLS BETWEEN 0 AND 2PI 
% 10/10/96 
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xX = direc; 
if x >= 2* pi 
x=x-2* pi; 
end; 
ix < 0; 
X=x+2*p1; 
end; 
out = x; 


LLZECEF.M 


function ecef=1 12ecef(x) 
% converts from geodetic latitude, longitude, elevation to % earth-centered earth-fixed Cartesian 
coordinates 
phi = x(1); lambda = x(2);h = x(3); 
% define semi-major and semi-minor earth axes a=6378137;b=6356000; 
% define auxiliary quantities f,e, and N 
f=(a-b)/a; e = f*(2-f); N = a/sqrt(1-(e*sin(phi))-2); 
% convert to Cartesian 
ecef = [(N + h)*cos(phi) * cos(lambda) 
(N + h)*cos(phi) * sin(lambda) 
(N*(I-e%2)+h) * sin(phi)]; 
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APPENDIX B 


The following Pages include manufacturer data-sheets from the DEL NORTE 


Technology, Inc. 


DEL NORTE Technology, Inc. e-mail: dnti@delnorte.com 


1) FLYING FLAGMAN® Global Positioning System Measuring Unit 


FEATURES/BENEFITS: 


12-Channel GPS Receiver 

Flow Rate/Temp/Humidity Monitoring 
10 Home Points 

100 Field Memory 

Polygon Field Entry 

Racetrack, Parallel Line and Squeeze to Middle Functions 
10 Position/Sec Updates 

Ground Speed in MPH, m/sec or knots 
Serial Port for control/logging 

Break Point Memory 

Resume to Point 

Automatic Reverse Direction Sensing 
Skipped Line Detection 

Selection of feet/meters 

Built-in Clock 

Data/Time Tag 

Remote Functions on Flight Stick 
On-Ground Self-test 


14] 


The new DEL NORTE Technology GPS Measuring Unit (GMU) provides a 
rugged, waterproof, and versatile GPS Receiver with multiple function capabilities and a 
range of options Set up is accomplished by using a menu driven Control and Display Unit 
(CDU). The Del Norte GMU can serve as either a Reference or Mobile GPS Unit. As a 
Reference Unit, the GMU, with its Establish Position capability, 1s able to derive site 
location to within one meter in a stand-alone mode over twenty-four hours and even better 
if left for forty-eight hours. 

After setup the CDU may be disconnected to leave the single rugged, waterproof 
unit on site. By adding a 'second processor’ option to the GMU, the Reference Station can 
also act as a differential monitor. As a mobile, the GMU may be interfaced to any of the 
user's data processing systems. For stand-alone operations the addition of the VGA board 
and optional second processor allow the GMU to drive a color display and incorporate 
third party software. 


SPECIFICATIONS: 


PERFORMANCE: 
e Differential Dynamic: 1 - 3 meters, 1 meter CEP 
e Update Rate : 0.1 seconds 
SYSTEM ENVIRONMENT: 
e Temperature: 0 - 55 degrees Celsius (32 - 130 degrees Fahrenheit) 
e Vibration: Tested to MIL-STD-810-C 
PHYSICAL/POWER CHARACTERISTICS: 


e System Weight: 24 lbs (10.88 Kg) e Super Trac Light Bar 

e CDU Dimensions e W 15.0 in (38.1 cm) 
e W 4.2 in (10.6 cm) e H2.01n (5.1 cm) 
e H/7.7 in (19.5 cm) e D5.5 in (14.0 cm) 
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eo) Diles in (Gecm) e Power Input : 11 - 32 VDC, 30 Watts 
e GMU Dimensions 

e W 8 in (20.3 cm) 

e H6.8 in (17.3 cm) 

e D 15 in (38.1 cm) 


SYSTEM SUPPLIED WITH: 
e GMU Mainframe e Differential Rx + Antenna 
e Control/Display Unit e Cables 
e Super Trac Light Bar e Shock Mounts 
e Remote Step Switch e Installation and Operation Manuals 
e GPS Antenna 
OPTIONS: 


e Differential Data Sources 

e Local Reference Tower 

e Portable Reference Station 

e FM Side-Band 

e Government MF Transmissions 

e Satellite Data Links 
e Data Logger and Post-Processing/Plotter Software 
e Map Trac Moving Map Display 


SUPER TRAC LIGHT BAR: 
e Pilot Programmable Displays e Left/Right Guidance 
e Digital HSI e Visible in Direct Sunlight 
e Line/Pass Number e Dimmer Control for Low Light/Night Flying 
e Distance Off-Track e Programmable Sensitivity 
e Distance-to-go e Course Deviation Indicators 


e Course Deviation e Programmable to nearest 0.1 of a degree 
e Area Sprayed 
e Custom Settings 


e Ground Speed 


2) MODEL 2012 PORTABLE REFERENCE STATION 
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FEATURES/BENEFITS: 


e Dedicated GPS Reference Station e 12-Channel GPS Receiver and Radio/Modem 
e Lightweight & Portable which can be in one box 
e Low-Cost e 'Self-Locating' 


The Del Norte 2012 is the low-cost portable reference station solution. To 
compliment the 2012, a small, lightweight, and inexpensive UHF Mobile Data Link 
Transmitter and Receiver is available to feed Differential GPS Corrections to your mobile 
GPS receiver. Del Norte has trained engineers and technicians available for field 
deployment anywhere in the world to provide assistance with first-time installations or 


emergency needs 


SPECIFICATIONS: 
PHYSICAL & POWER CHARACTERISTICS: 
e Weight 5 Kg e Operating Temperature 
So Dacncion< 0 - 50 Degrees Celsius 
W 30 cm e Power 
9 cm 11 - 14.5 vDC 
D 34 cm 
SUPPLIED WITH: 
e 12-Channel L1 Receiver e Necessary Cables and Connectors 
e UHF Transmitter and Modem e Operator's Manual 
e Patch GPS Antenna w/Ground Plane e Comprehensive Product and Firmware 
e Broad Band Omni Antenna Maintenance for One Year 


e Training Available 
SYSTEM OPTIONS: 
e 215UL Mobile Data Link Receiver 
RADIO SPECIFICATIONS: 


e Frequency 
458.xxx MHz, 0.5 Watt (UK deregulated), options include other UHF, VHF, and 
HF frequencies on request 


e Transfer Times 
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W/RTCM Format 
Typical w/4 satellites, 0.4 seconds 
Typical w/12 satellites, 1.0 seconds 


Baud Rates 
300 to 9600 
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